Spin nematic state as a candidate of the hidden order phase of URu2Si2 
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Motivated by the recent discovery of broken four-fold symmetry in the hidden order phase of 
URu2Si2[R. Okazaki et al., Science 331, 439 (2011)], we examine a scenario of a spin nematic 
state as a possible candidate of the hidden order phase. We demonstrate that the scenario nat- 
urally explains most of experimental observations, and furthermore, reproduces successfully the 
temperature dependence of the spin anisotropy detected by the above-mentioned experiment in a 
semi-quantitative way. This result provides strong evidence for the realization of the spin nematic 
order. 

PACS numbers: 



The heavy fermion compound URu2Si2 exhibits a sec- 
ond order phase transition at Tjjo ~ 17.5 K. In spite of 
long-standing enormous efforts in experimental and the- 
oretical studies [l- 14 1, the order parameter of this phase 
transition has not yet been identified. The enigmatic 
features of this so-called "hidden order(HO)" phase are 
described as follows: (i) Despite large anomaly in ther- 
modynamic quantities and drastic reconstruction of the 
Fermi surfaces at T = Tho, there is neither conventional 
magnetic order nor the change of the crystal structurefl- 
Isl. [isl - flij . (ii) However, under applied pressure, an an- 
tiferromagnetic (AF) state with large moment appears, 
and more surprisingly, the Fermi surfaces in the AF or- 
dered state are almost the same as those found in the HO 
phaseliiil. 

Recently, an experimental breakthrough for this issue 
was achieved by Okazaki et aH25], who found sponta- 
neous symmetry breaking in the spin space at T < Tho- 
They reported that the anisotropy of the spin suscepti- 
bility in the xy-plane, which is measured by the quan- 
tity Xxy = {SxSy), becomes nonzero below Tho- Since 
URu2Si2 is tetragonal with four-fold symmetry at T > 
Tho 7 and the phase transition at T = Tho does not ac- 
company any lattice distortion, it is reasonable to ex- 
pect that this symmetry breaking is an essential feature 
of the HO phase, which imposes a crucial constraint on 
possible candidates of the HO parameter. Motivated by 
this experimental observation, in this letter, we discuss 
a possibility that a spin nematic (SN) state is realized 
as the hidden order in URu2Si2. The SN phase is a 
state with circulating spin currents, but with no magnetic 
moment 26l-l28l| . The circulating spin currents break spin 
rotational symmetry, leading to spin anisotropy, with- 
out breaking time-reversal symmetry. We demonstrate 
that the above-mentioned features of the HO phase are 
naturally understood within the scenario of the SN or- 
der, and furthermore, that the temperature dependence 
of the spin anisotropy spontaneously generated in the SN 
state successfully explains the above experimental obser- 
vation in a semi-quantitative way, providing strong evi- 
dence for the realization of the SN state as the HO phase 



of URu2Si2. 

We, first, present a mean field analysis for basic prop- 
erties of the SN state applied to the case of URu2Si2. 
The SN state is induced by nesting of the Fermi sur- 
face as discussed in refs.^26, 27]. In fact, the recent band 
calculations for URu2Si2 based on an itinerant /-electron 
picture revealed that there are one electron band denoted 
as Efej^ and one hole band denoted as Efcj [H, 29], which are 
nested to each other via the nesting vector Qo — (0, 0, tt): 



i.e. Ek+Qoi = —£k2 |30|,|3l[. It is noted that Qo is equiva- 



lent to the ordering vector of the large-moment AF state 
which is observed under applied pressure [20|. Moreover, 
it was pointed out that several experimental results sug- 
gest that the energy gap opens on the Fermi surface below 
Tho Hj 32-3^, which is consistent with the gap genera- 
tion due to the nesting of the Fermi surface. Because of 
these reasons, we employ the scenario that itinerant /- 
electrons undergo the transition to the SN state triggered 
by the Fermi surface nesting. We will discuss the micro- 
scopic origin of this instability later. Since the ordering 
vectors for the large-moment AF phase and the SN phase 
are the same, the feature (ii) mentioned above is natu- 
rally understood within our scenario. The SN phase is a 
spin-triplet electron- hole pairing state 2^ 53] , and hence 
the order parameter of the SN state for the effective two 
band model is 



(4l<TCfc+Qo2<T') = di2{k) ■ CTa 



(1) 



where c]^^^^ [ckaa ) is a creation (an annihilation) operator 
for an electron in the band a = 1, 2 with momentum fc, 
spin (T. di2{k) is a vector, the direction of which is paral- 
lel to the spin quantization axis of the SN order. Because 
of time-reversal invariance in the SN state, we have the 
condition <il2{^k) = —di2{k). Furthermore, we impose 
inversion symmetry, since there is no indication of broken 
inversion symmetry in URu2Si2 from experiments. Then, 
it follows that dl2{k) = —di2{k). Also, we assume that 
dab{k) is symmetric with respect to the exchange of the 
band indices a and b. This assumption will be plausi- 
bly justified by a microscopic argument given later. The 
mean field Hamiltonian for the SN state of the effective 
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^iiF "with the kinetic 



two-band model is Hmf = T~Lmf 
energy term r6y^^ = T.k,aY.a=i.,2^'^'^'^kaaCkaa, and the 
particle-hole pairing term, 



Ckba 



(2) 



k,<7,a' 



For URu2Si2, the crystal structure of which has the 
Dih symmetry, the anisotropy of the spin susceptibil- 
ity which breaks four-fold symmetry down to two-fold 
symmetry in the xy-plane implies that the order param- 
eter belongs to two-dimensional (2D) representation of 
the D4fi symmetry, and in the hidden order phase, four- 
fold symmetry in the 2D space is spontaneously bro- 
ken. Then, from the symmetry properties of <ii2(fc) dis- 
cussed above, a possible candidate is the Eg state with 
di2{k) = i{Aikyk^,A2k^k^,0) (or i{Aikj;k^, A2kyk^,0)) 
for small |fc|. Here, the real parameters Ai and A2 are 
determined from the self-consistent gap equation. 

According to the experiment^], the axis of the Ising- 
like anisotropy, which is spontaneously generated in the 
HO phase, is parallel to the (l,l,0)-direction. Since our 
toy model has continuous rotational symmetry in the xy- 
plane, the analysis given here is applicable to URu2Si2 by 
rotating the principle axes by 7r/4 around the z-axis, i.e. 
x' — -^{x~y) and y' — -^{x + y). Following the exper- 
imental observation, we assume that four-fold symmetry 
in the x'y'-plane is spontaneously broken, resulting in the 
state with di2{k) \\ (1, 0, 0) or (0, 1, 0) in this rotated spin 
frame. The direction of di2{k) is determined by the de- 
tail of the electronic structure and spin-orbit interaction. 
However, most of the following results do not depend on 
it. For the tight-binding model, we choose 

di2{k) = {iAi(f>{k),0,0), (j){k) — sink fj,' sin k^, (3) 

where /i' = x' or y' . 

At this stage, we note that the electron-hole pairing 
term ([2]) is nonzero only when the nesting vector Qq and 
the momentum dependence of dab{k) fulfill the following 
relation; e~''5°''' — e*'^°''^ ^ for r^, rj satisfying Aij = 
J2k(t>ik)e-'''^'-'-'-^'> ^ 0. For © and Qo = (0,0, tt), 
this is actually fulfilled. It is instructive to compare this 
property of the SN phase with the unconventional spin 
density wave (USDW) state considered in refs. 0,0, the 
order parameter of which is also given by Eq.([T]) but with 
di2 a real even function of fe, because of broken time- 
reversal symmetry. It is easy to see that, if Qo = (0, 0, tt), 
and the momentum dependence of di2(fc) is the same as 
(|3]), the particle-hole paring term ([2]) vanishes, and hence, 
this type of the USDW can not be realized. It is noted 
that the USDW is more stabilized than the SN phase, if 
the order parameter dab{k) (and hence, the particle- hole 
pairing interaction) is anti-symmetric with respect to the 
exchange of the band indices a, h, or if higher harmonics 
of the order parameter such as 0(fe) = sin2fc^/ sin2kz is 
allowed. 
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FIG. 1: (a) Solid line: anisotropy of the spin suscepti- 
bility Ax = Xx'x' — Xy'y' in the SN phase versus T/Tuo- 
The magnitude of Ax is normalized by the value at T = 0. 
Circle: experimental data quoted from [2^. Ai(0)/Tho = 
2.6. (b) Imx^^(Qo,'^) versus lj/Ai. ^b-H"z/Ai = O(solid), 
O.l(dotted), 0.2(dashed), 0.3(dot-dashed), 0.4 (double- dot- 
dashed), (c) Imx^^(Qi,i^) versus cj/Ai. /iB-ffz/Ai = 
O(solid), O.l(dotted), 0.2(dashed), 0.3(dot-dashed). 



We, now, calculate the anisotropy of the uni- 
form spin susceptibility, which characterizes the SN 
order, and was experimentally detected in ref . (25| . 
The anisotropy in the x'y'-plane for our model 
is obtained as, AxiT) = Xv'y'{T) - Xx'x'{T) - 



tanh - 
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where Ef. 



2Tcosh2 

^, and the momentum sum EL is taken 



V4i + |di2(fc)l 

over the nested part of the Fermi surface satisfying the 
condition Sk+Qai — —£k2- Also, we assume the BCS 
mean-field-like T-dependence of the amplitude of the gap 
Ai (T). To make a semi-quantitative comparison with the 
experimental data without referring to the details of the 
band structure, we consider a ratio Ax(T)/Ax(0), for 
which it is expected that effects of specific band struc- 
tures approximately cancel out. There are two parame- 
ters in this calculation; one is the ratio of the energy gap 
at T = to the transition temperature, i.e. Ai(0)/Tho, 
and the other is an overall normalization factor of the 
magnitude. According to the recent scanning tunneling 
microscopy (STM) measurement, the magnitude of the 
gap opened on the Fermi surface via the hidden order 
transition is Ai(0) ~ 4 meVjH, HI. Since Tro = 17.5 
K, we choose the parameter as Ai(0)/Tho — 2.6. Then, 
there is only one fitting parameter, i.e. the overall nor- 
malization factor. We choose the normalization factor 
to fit the theoretical result at T = 6 K to the experi- 



3 



mental data at the same temperature. The calculated 
result of this one-parameter fitting is shown in FIG[lJa). 
In spite of simplicity of the model, the theoretical result 
is surprisingly in good agreement with the experimen- 
tal observations, at least above T/Tho ~ 0.15. A slight 
discrepancy at low temperatures may be ascribed to the 
superconducting transition. This result provides strong 
evidence for the realization of the SN state as the HO in 
URu2Si2. In FIGlUJa), the T-linear behavior of Ax(T) 
at low temperatures is seen, which is raised by the exis- 
tence of the line node of the order parameter ([3]). The 
experimental data is also consistent with this behavior 
above T/Tuo ^ 0.15. 

The above scenario also has important implications for 
magnetic properties of the HO phase. According to neu- 
tron scattering measurements, there is a longitudinal spin 
fluctuation with the wave number q = Qq, which exhibits 
an excitation gap in the HO phase [sJl . Remarkably, the 
magnitude of the spin gap ^ 1.6 meV is much smaller 
than the single particle energy gap Ai ~ 4 meV ob- 
served via the STM measurements 32,33|. Furthermore, 
the spin gap increases notably with increasing a magnetic 
field applied along the z-axis, Hzi_e.g. the spin gap for 



Hz = 17 T reaches to 2.5 meV [SJ]. These properties 
are well explained by the present model. Using the mean 
field Hamiltonian 'Hmf and the random phase approx- 
imation, we calculate the longitudinal spin correlation 
function x'^Hq,^) = dt{[S''iq,t), S''{-q,0)])e"^* 

for q = Qo, which is dominated by the transition be- 
tween the band Ski and the band ek2 [IB]. In FIG.[lJb), 
the imaginary part of x^^iQoj is plotted as a function 
of the frequency w/Ai, which indicates the spin excita- 
tion gap smaller than the single-electron gap Ai. As the 
magnetic field increases, the excitation gap increases 
substantially. These behaviors are understood as a result 
of the existence of line nodes of the gap |<ii2(fe)|. The line 
nodes allow low-energy spin excitations to develop below 
Ai. When the system is in the vicinity of magnetic crit- 
icality, a gap structure appears for w < Ai, which is also 
sensitive to applied magnetic fields, as shown in FIG. 
[TJb). These results are in good agreement with the ex- 
perimental observation obtained in ref.js^l- Conversely, 
these neutron scattering data strongly imply the exis- 
tence of line-nodes in the energy gap. On the other hand, 
these experiments also revealed that in addition to Qo, 
there is a longitudinal incommensurate spin fiuctuations 
with Qi = (1.47r,0, 0) ^1], which may be attributed to 
the nesting between the hole band ek2 and another hole 
band denoted as Sks, as suggested from the band calcu- 
lations ^,130f . According to the recent experiments, the 
magnetic excitation for the incommensurate Qi has an 
energy gap ^ 4 meV below Tho UM, ISJ, l37| . In our sce- 
nario, the gap ~ 4 meV opens in the hole band ek2 in the 
SN state, while it does not in the hole band £fc3. Thus, 
the magnetic excitations due to the transition between 
these two hole bands should exhibit the excitation energy 



gap ^ 4 meV at the wave vector Qi rather than 2x4 
meV. This explains the above-mentioned experimental 
result. To demonstrate this, we calculate the spin corre- 
lation function x^^(<5i, ui) for the mean field Hamiltonian 
(36| . In this calculation, we postulate that hot spots on 
the hole band 2 for this spin fluctuation are located away 
from the gap-node of |di2(A;)|; i.e. |(/i2(A;)| ^ for fc 
satisfying ek+Qi2 ~ —£k3- Then, the sharp gap edge ap- 
pears at w = Ai in Imx^^(Qi, i^), as shown in FIG.[l{c). 
In the case with magnetic fields Hz, the gap edge shifts 
to w « Ai -I- 2iJ2/Ai Thus, the field dependence of 
Imx^^(Qi,a;) is much weaker than that oiIm.x^^{Qo,uj) 
(FIG.[lJc)). This is also qualitatively in agreement with 
the experimental results [sj, [s^]. 

Some remarks are in order: (i) In the SN state, there 
are staggered circulating spin currents [H, [23|. When 
a uniform magnetic field parallel to (1, 1,0) or (1,-1,0) 
in the original frame is applied, and electron spins are 
polarized, the spin currents lead to the staggered orbital 
currents, which induce staggered moment, as in the case 
of the orbital current state 0. In NMR experiments, 
broadening of spectra under applied fields was observed 
Til 17 1 . Its origin may be attributed to the staggered 



circulating currents of the SN state. The magnitude of 
the induced staggered field raised by the circulating cur- 
rents is estimated as ^ 1 Oe for the applied field 7J ~ 4 



T, which is consistent with the experimental results|17|. 

(ii) Note that the USDW state with the spin quanti- 
zation axis parallel to (1,±1,0) also exhibits the same 
spin anisotropy as the SN state, and may explain the ex- 
perimental data of ref.ji^. At this stage, we can not 
thoroughly exclude the possibility of the USDW as a 
candidate of the HO, though, as mentioned before, the 
particle-hole pairing interaction anti-symmetric with re- 
spect to the band indices, the origin of which is micro- 
scopically unclear, is required for its realization. 

(iii) The SN order parameter couples to a magnetic 
field H a.s ^ (di2 • H)^ in the free energy. This im- 
plies that the nonlinear susceptibility xs for a magnetic 
field H satisfying di2 ■ H exhibits a discontinuous 
jump at T = Tho for the SN order ([3]). However, such 
an anomaly was experimentally observed not for the in- 
plane field, but for || z-axis |3j. The discontinuous 
jump for H \\ z-axis may be explained by postulating 
that the SN vector di2 is not confined in the xy plane, 
but rather has a nonzero out-of-plane component. This 
scenario may also explain the change of the slope of the 
linear spin susceptibility at T = Tho for H || z-axis. 
However, the absence of the jump for the in-plane field 
in the experiment remains to be resolved. It is desirable 
to re-examine the measurement of the nonlinear suscep- 
tibility using recent high-quality samples. 

(iv) According to specific heat measurements in mag- 
netic fields parallel to the z-axis H^, the specific heat 
jump AC at T = Tho is almost constant as a function 
of Hz, though the transition temperature decreases as 
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Tno{H^) - rHo(O) oc -Hi [Hill. These behaviors are 
explained by the SN scenario. The decrease of Tho is un- 
derstood as a result of the decrease of the particle-hole 
pairing interaction due to the magnetic field [s^. Also, 
the absence of the change of AC under applied fields is 
explained by taking into account the selfenergy correc- 
tion arising from interactions with AF spin fluctuations 
[sif . Within this scenario, the specific heat jump is given 
by AC ^ ^^ThOj where ^ is the correlation length for AF 
spin fluctuations. If the system is in the vicinity of the 
AF critical point, i.e. ~ AC is constant, though 
Tho is decreased by the magnetic field. 

(v) Experimental studies reported that under applied 
pressure, the discontinuous phase transition from the HO 
state to the large moment AF state occurs [i^, while 
the Fermi surfaces in both phases are almost the same 
[3 21, 23, 24 1 . According to our scenario, there are 
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competing AF spin fluctuations and fluctuations toward 
the SN phase transition in the system, (see below and 
Supplemental Material for details) Thus, applied pres- 
sure induces the change of effective couplings between 
electrons and these fluctuations, resulting in the transi- 
tion to the AF phase. Since the order parameters of these 
phases have the different symmetries that are not com- 
patible with continuous phase transition, the type of the 
phase transition is first-order, which is in agreement with 
experiments 17, H^]. We stress that the reconstructed 



Fermi surfaces in the SN phase and the AF phase are the 
same in our scenario, because both of them are recon- 
structed by the same nesting vector Qq. This is also in 
agreement with experimental observations 3 21-24|. 

Finally, we discuss microscopic mechanism of the SN 
order for URu2Si2 [s^]. We consider a scenario that or- 
bital fluctuations associated with the Fermi surface nest- 
ing yields the SN order. The minimal model consists of 
the three bands Ska with a = 1, 2, 3, and mutual Coulomb 
interaction between electrons. The SN order considered 
here arises from the particle-hole pairings in the d-wave 
channel which are formed between between the electron 
band 1 and the hole band 2. The Fermi surface nesting 
with the nesting vector Qo between these bands leads 
to this instability. An effective pairing interaction in 
the (i^i'z-wave channel is mediated via orbital fluctuations 
arising from the Fermi surface nesting between £^3(1) and 
£k2 with the nesting vector Qi (Qo)- If the orbital fluc- 
tuations are sufficiently strong, the SN order is stabilized 

In conclusion, we have demonstrated that the SN state 
induced by the Fermi surface nesting is a promising can- 
didate of the HO phase of URu2Si2, since it successfully 
explains most of experimental observations including the 
recent experiment on four-fold symmetry-breaking. 
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in these bands. The Hamiltonian is 

+ 5^ ^criCr2CT30-4''l+qa(Ti'^l'-qf,(T2'"'=''='^3Cfcd(T4,(4) 

k,ki' ,q 

with a, b,c,d— 1, 2, 3. The SN order occurs for electrons 
in the band 1 and the band 2, i.e. |<ii2(fc)| = |<i2i(fc)| ^ 0, 
and \dab{k)\ = for other a, b. 

Following the results of the recent band calculation [sij, 
we postulate that the three bands, Cfci, ek2, and Sks, have 
the orbital character with = ±5/2, ±3/2, and ±1/2, 
respectively, and that each of the single electron states in 
the a-th band with spin a denoted as \a a) approximately 
corresponds to the orbital state \Jz)- That is, 

|1 1) = I J. - 5/2), lU) = 1-5/2), (5) 
|2t) = |3/2), 12;) = 1-3/2), (6) 
|3t) = |l/2), |3;) = 1-1/2). (7) 

The interaction term of (|H) conserve the total Jz in each 
scattering process. We, furthermore, assume that in 
the interaction term of interactions between elec- 
trons with anti-parallel spins dominates over those be- 
tween electrons with parallel spins, and neglect the lat- 
ter. Then, nonzero coupling constants for the interaction 
between electrons in the band 1 and the band 2 are, 



T raabb 



T raabb 



jjabba 



Jjabba 



(8) 



with (a, 6) = (1,2) or (2,1). For simplicity, we assume 
that all of these coupling constants are the same, and 
equal to U\2- Since we are concerned with effects of an 
applied magnetic field H , we also add the Zeeman term 
—giJ^BS^H to the Hamiltonian Assuming the per- 
fect nesting condition Sk+Qai = ~£k2, we write down 
the non-interacting single-electron Green function for the 
band 1 with spin a in the SN state, Gkiai^n), in the fol- 
lowing forms. 



SUPPLEMENTAL MATERIAL 



A. Spin correlation functions in the spin nematic 
state 



GklaiSn) 



■ Efcl 



Ha 



(9) 



where h = g^iBH and Eka- = \/(efci - s^h)^ + \di2{k)\^ 
with So- — +1 for cr =t and —1 for a The Green 
functions for the band 2 are 



In this section, we present the calculations of the dy- 
namical longitudinal spin correlation functions x'^^iQi'^) 
in the SN state, particularly focusing on the modes with 
q = Qo or Qi, for which spin fluctuations are promi- 
nent. The minimal model consists of the three bands, 
£fei, £fc2, and £^3, with the particle- hole pairing term, 
and the mutual Coulomb interaction between electrons 



Gfe2t(£n) = — Gfel^( — £„), 
Gfc2i(£n) = —Gfcl-f (-£«)■ 



(10) 
(11) 



The non-interacting anomalous Green functions gener- 
ated by the SN order are given by. 



Jo [^Sn) 



2~rE^ ' 
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i\di 2{k)\ 



The longitudinal spin correlation function enhanced by 
nesting between the band 1 and the band 2 is 



not open in the band 3 in our model, we obtain within 
the random phase approximation, 



(0)/ > 



1 - U23.X2i{(l^^n) 



(20) 



where 



= X?22l(g:^n) +X2112(g:^^«) 
+ Xm2(g;^«) +X212l(g:'^") 
= 2[xl22l(9,W„) +X2?2l(9,W„)], (14) 



dT{T{S:,{q, T)S:,{-q, 0)})e*-"-(15) 



where 



X23 



E 

k,(T 



tanh ^ + tanh ''"i+f" 



Eka + Ski + Sah 



■w+{k) 



tanh % - tanh ^ ^ 
^ ^^——w^{k) 

lUJn - Eka + £kl + Sah 



with ^ 5Efe(4+qaTCfebt - 4+qa;Cfcb;). In the 

second equality of (HH), we used the relations xf|2i = 

X2II2 S-nd X12I2 = X2121- 

Within the random phase approximation, xf | is given 

by 



with 



(21) 



(22) 



(0)/ ^ 
X12 (9>^") 



1 TT (0)/ 

1 - 1/12x12(9.'^" 



(16) 



Imx^^(Qi,w) is approximated by lmx23^ {Q 1 , t^) ■ 

Plots of Imxi2''^(Qo, i^) and lmx23^ {Q 1 1 func- 
tions uj are shown in FIG. 1(b) and (c), respectively, in 
the main text. 



where xi2''(9''^n) the spin correlation function in the 
case without Coulomb interaction. For q — Qq, it is 
given by, 

T 

X12 (Qc^n) ^^EE/t E ^fclo'(^"j)G'fc2(T(em + W„) 

m k <7=^;,i 



tanh ■ 



E 



tanh 



(i/fc^ + Ek^Y - iliOnY 



^ ~ tanh ^ 



2T 



2T 



(i;fct - EkiY 

where 



{iUn)' 



■{Eki - Eki)u^ik) 



(17) 



B. Microscopic mechanism of the spin nematic state 

In this section, we discuss microscopic mechanism of 
the SN order for URu2Si2. We consider a scenario that 
spin fluctuations and orbital fluctuations associated with 
the Fermi surface nesting with the nesting vectors Qq 
and Qi give rise to the phase transition to the SN order. 
Our argument is based on the three band model ([4]). 

The spin fluctuations and the orbital fluctuations with 
the propagating wave numbers Qo and Qi mediate the 
effective interactions between electrons. 



u±{k) 



1 ± 



'-kl 



\d,2{k)\'-h^ 



Ek-fEki 



E 



O"! fT2 Cr3 (74 

C/cd(T4(l23) 



(18) 



Imx^^ (Qo I <^) is approximately given by the retarded spin 
correlation function lv[iXi2^'{Qo,^)- 

In a similar manner, we can obtain x^^(Qi,w), which 
is dominated by the longitudinal spin correlation function 
enhanced by nesting between the band 2 and the band 3, 

X2f(g,^n) - 2[x^|32(q,W„) +X2l23(9>'^n)]- (19) 

From Eqs.® and ([7]) and the conservation of Jz in inter- 
action processes, we see that nonzero coupling constants 
of the interaction between electrons in the band 2 and 
the band 3 are given by Eq.([8]) with (a, 6) = (2,3) or 
(3,2). We assume that all of them are equal to U2j,- To 
simplify the calculation, we also impose the perfect nest- 
ing condition ek+Qi2 = — £fc3- Since the energy gap does 



with q = (q, u). 

From Eqs.®, ([5]), and (O, and the conservation of the 
total Jz in each scattering process, the effective interac- 
tion that mediates the particle-hole pairing characteriz- 
ing the SN order is approximately given by. 



p2211 



(9) 



^^^i;("H'/)-ffj#?ff|?II[xj?(g) + xI!(<z)] 
-5nn5!tn[xn(9) + xu(9)]- (24) 



Here, the second and third terms are the interactions me- 
diated via the orbital fluctuations which is characterized 
by the orbital susceptibility. 



xt{<i) 



dT{T{NabAQ, r)Nbaa{-q, 0)})e'^"\ (25) 



7 



and 



where Naba{q) — Y^k '^k+qaa^kba ■ ^tttt' ^Ittl' ^TUT 

^l^ll are effective coupling constants for interactions be- 
tween the orbital fluctuations and electrons, which pre- 
serve Jz- The first term of (|24|) is the other interac- 
tion processes. The Fermi surface nesting between the 
band 1(3) and the band 2 with the nesting vector Qo 
(Qi) enhances the orbtial- fluctuation- mediated part of 
((24)) . while the first term of (|24p is not much affected 
by it. In Ea. (l24l) . we have neglected the interaction pro- 
cesses that are mediated by orbital fluctuations associ- 
ated with x"J(9) and x"^j because these orbital fluctu- 
ations accompany transverse spin fluctuations, and, ac- 
cording to neutron scattering measurements, transverse 



spin fluctuations are fully suppressed in URu9Si9.ll 5 1 
The second and third terms of (f24| give rise to the 
particle-hole pairing interaction in the d^'z-wave channel 
Vd ^ nil' '/'(fc)0(fc')r?JI|K,fc - fe'). Here EL'- de- 
notes the momentum sum over the Fermi surfaces. Since 
(j){k)(t){k + Qrn) < with to = 0, 1, and the orbital 
susceptibility (|25|) is significantly enhanced at q = Qi (or 
Qo), the sign of Vd is positive, which leads to the d-wave 
particle-hole pairing characterizing the SN state. 



In addition to (f24)) . there is another important inter- 
action channel which is related to the instability toward 
antiferromagnetism with staggered magnetization paral- 
lel to the 2;-axis. This instability is also driven by the 
Fermi surface nesting between the band 1 and the band 
2, as in the case of the SN order j30|. We denote this ef- 
fective interaction as T^J!f.^{q). Note that r^^|j_((7) is not 
enhanced by the orbital fluctuations associated with the 
orbital correlation function ([^5)) . Thus, if the orbital fluc- 
tuations are sufficiently strong, the effective interaction 
((24| gives rise to the instability toward the SN state. The 
transition temperature of the SN order within a mean 
field approximation is determined by the self-consistent 
gap equation. 



k 



cf'ik). (26) 



When there is an applied magnetic field parallel to the 
2;-axis, the orbital susceptibility p5|) is reduced. For a 
small field i?^, the nonzero lowest order correction is 
'^Xctct(9) —co((7)i/f with co(g) a function of q. The mag- 
netic field parallel to the z-axis affects the self-consistent 
gap equation p6)) only through the effective pairing in- 
teraction. Thus, the transition temperature in the case 
with Hz behaves like Tho(^^z) ~ ?Ho(0) — a^Hl with oq 
a positive constant. This field dependence of the tran- 
sition temperature is actually observed in experiments 

fii. 



C. Specific heat jump 

In this section, we discuss the specific heat jump at the 
HO transition point. According to specific heat measure- 
ments under magnetic fields applied along the z-axis, the 
specific heat jump AC at T = Tfjo is almost constant as 
a function of the magnetic fields, though the transition 
temperature decreases as Tho(^^z) ^ T'Hn(O) — aoH^, as 
mentioned in the previous section [bsI . l39l |. The constant 
specific heat jump is unusual for the itinerant f-electron 
system, since in the Fermi liquid state, the entropy gen- 
erally decreases as T decreases. However, here, we show 
that this remarkable feature is explained by taking into 
account the selfenergy corrections arising from AF spin 
fluctuations. We start from the Ginzburg-Landau (GL) 
free energy in the vicinity of the HO transition point, 

= aip'^ + btp"^, with tp the order parameter for the SN 
state. We have chosen the spin quantization axis parallel 
to di2. Thus, Tp — Ai is a scalar order parameter. The 
speciflc heat jump at T = Tho is given by 



AC: 



HO 
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da 
df 



(27) 



T=Tho 



The coefficients of the GL free energy a and b are calcu- 
lated from the normal Green functions. 



Gfci(j(e,i) — 



1 



71 ^kl 



with 



Zki^n) = 1 — -. Sfc(e„), 



(28) 



(29) 



and S]fc(e„) the selfenergy, and Gk2ai£n) which is ob- 
tained from the relations ([TOl) and (fTTj) . The results are 



a=--n{Qo) 
9 



(30) 



with 



n(Qo)--T J2 I]</''(fc)Gfcia(e„)Gfe+Q„2-.(e„)(31) 
and 

b = T J2 E'^'('^)[^'^i-(^«)^'^+Qo2-^(en)]'- (32) 

|n|<nc k 

Here ~ E^/T with Ec an energy cutoff for the effective 
particle-hole pairing interaction g. For simplicity, we ne- 
glect fe-dependence of Zfc(e„). Then, from egs. ipT]) . (PT|) . 
and (|32]). we obtain. 



AC«7V(0)Z(^,)Tho, 



(33) 



where N{0) is the density of states at the Fermi level. 
The selfenergy corrections Z{Ec) mainly stem from in- 
teractions with AF spin fluctuations and orbital fluctu- 
ations. Note that the AF spin fluctuations expressed by 
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(IT4| and (fT9|) and the orbital fluctuations characterized 
by ((25)) are not independent, but deeply correlated. The 
former contributions generally enhance the latter, and 
vice versa. We postulate that the normal selfenergy cor- 
rections are mainly due to the interaction with the AF 
spin fluctuations with Qq, while the particle- hole pairing 
interaction (l24l) is dominated by the fluctuations with 
q ~ Qi; i.e. in Eg. ([24| . the coupling strength of the 
second term is much smaller than that of the third one, 
though the correlation length for the spin fluctuations 
with Qq is larger than that with Qi. We also assume 
that the characteristic energy scale of the AF spin fluctu- 
ations with q = Qq is sufficiently larger than the Zeeman 
energy of the applied fields, and the AF spin fluctuations 
with Qq is not much affected by the magnetic fields, at 
least, in the vicinity of T ~ Iho- The most important 
effect on Z{sn) in the vicinity of the transition tempera- 
ture is the quasiparticle damping 7. If the system is suf- 



ficiently clean, the quasiparticle damping raised by the 
interaction with the three-dimensional AF spin fiuctua- 
tions with Qq is 7 ^ ^2 [si + (7rT)2] with ^ the correlation 
length for the AF spin fiuctuations [i^ . Then, from (|33|) , 
we have AC ~ C^Iho up to a constant factor. If the sys- 
tem is in the vicinity of the AF quantum critical point, 
i.e. ~ l/T, AC does not depend on Tho- This result 
provides a possible explanation for the above-mentioned 
experimental observation that the specific heat jump at 
T = Tho is not changed by the magnetic field, though 
the transition temperature is decreased by the magnetic 
field as TnoiH^) ~ Tho(O) -aoF^^ [H,!!^. We note that 
the above scenario of the constant specific heat jump is 
valid as long as Tho(^^z) is sufficiently large compared to 
the change of the characteristic energy of the spin ffuctu- 
ation raised by the magnetic field. Thus, when Tho(^^z) 
is decreased too much by the applied field, the deviation 
from the above result may appear. 



